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I. INTRODUCTION 

In response to a tensile stress glassy polymers can exhibit 
crazing which is a phenomenon of considerable practical and 
theoretical importance in the use of many plastics. Since the 
first important paper [1] dealing with the craze was published, 
a great de%l of knowledge has been accumulated on the development 
of crazing in glassy polymers in a tensile stress field [2-9] . 

For predicting craze growth in polycarbonate, Verheulpen-Heymans 
and Bauwens [7] have presented a two-dimensional model for craze 
growth. They obtained that the craze length varies linearly with 
respect to the log of time and hence craze growth velocity de- 
creases as time increases, in good agreement with their experi- 
mental result under constant uniaxial tensile load. Argon et al 
[8] have proposed new mechanisms for craze initiation and growth. 

In their model, the craze growth velocity is assumed to be pro- 
portional to the constant craze tip opening velocity. Then they 
conclude that the craze length increases linearly with time. They 
also get good agreement between theoretical and experimental results 
on PMMA and PS. 

It has been pointed out that both of the above models were 
somewhat similar in many ways and contained similar limitations. 
While they satisfy specific data but are incompatible and 
neither are capable of predicting the effect of craze interactions 
[ 10 ] . 

Recently, a time dependent theory of crazing behavior on 
polymers [11] based upon linear viscoelasticity theory and first 
law of thermodynamics for predicting the craze length as a function 
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of time, has found that a craze may propagate with increasing 
velocity for a thin sheet of viscoelastic medium which contains a 
single craze under uniaxial creep conditions. 

It is difficult to draw conclusions for one against the 
others, because there are several potential sources of discrepancy 
among the results of [71 , [8] and [11] such as tensile stress 
applied, material used and environment exposed. However, an 
important factor which may dominate craze growth behavior ia the 
interaction of the neighboring crazes* 

In this paper, an attempt is made to demonstrate how craze 
growth behavior is affected by the craze density (which is de- 
fined as the total number of crazes per unit surface area) 
surrounding that craze. The measured craze length as a function of 

time will be compared with some calculated results of PS at room 
temperature. 


II. EXPERIMENTAL 


As shown in Figure 1, the development of multi-crazes in 
a simple tension specimen is usually on the surfaces of the speci- 
men. In order to measure the changing length of the crazes, an 
experimental set up for determining the craze tip position and 
thus the craze length as a function of time is constructed. 

Sheet specimens of polystyrene with a neck region have been 
cut and secured onto the clamps of a creep testing machine for 
loading. As shown in Figure 2 while the specimen was being 
stressed, a microscope was used for viewing the surface conditions 
of the stressed specimen. Using a He-Ne laser the surface crazing 
was projected onto an image screen for observation and measurement. 
All the creep tests were conducted at a constant room temperature 
of 21°C under a constant applied load. Usually a short period 
of time elapsed before crazes occurred. Starting at this craze 
initiation time the craze lengths were monitored. The craze lengths 
were enlarged 400 times and measured to 0.1 mm at regular intervals 
with the aid of a travelling microscope. 

The eight inch-long specimens were cut from a 10-mil thick 
biaxially oriented polystyree sheet. The gage section was half 
inches wide and two inches long. Under load after inception of 
crazing , the number of crazes (per unit area) as well as the 
lengths of individual crazes were observed to increase. Figure 
3 shows the total number of craze per unit area as a function of 
time for a typical specimen under the load corresponding to a 
creep stress of 34.4 N/mm 2 at the gage section. Figure 4 
shows the increase in the individual craze lengths as time increases 


Biaxially oriented polystrene specimens were used because of 
the ease of developing the crazes. The density of the craze was 
found to be adequate for quantitative considerations. During the 
course of this investigation sufficiently large areas have been 
taken to insure the statistical representation of the crazing 
behavior. It was recognized that both the lengths and the number 
of crazes increase as a function of time. As shown in Figure 5, 
in a local region, out of many crazes, four were selected for 
identifying a common behavior. The craze lengths were plotted 
as a function of time. The similarity of their behavior is 
quite clear. More interesting is that by normalizing the craze 
length (current craze length divided by its initial length) all 
the individual curves reduce to a single one. This is shown in 
Figure 6, a single master curve characterizes the time dependent 
behavior of all the crazes developed in that region. In other 
words' this single master curve represents the typical behavior of 
a family of crazes in a local neighborhood under a specific loading 
situation as stated earlier. 
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III» REDUCTION OF THE MULTI -CRAZE PROBLEM 

Among all the parameters describing the crazing phenomenon, 
two important ones are the total number of crazes appearing on 
the surface of the stressed medium N(t) at time t and the length of 
each (ith) individual craze c i (t). The former is associated with 
the craze initiation while the latter is related to craze growth. 

Both parameters are highly dependent on the viscoelastic behavior 
of the media. 

In dealing with crazing, like any other physical or mechanics 
problems, continuum models are utilized. Aside from the regular 
field equations such as equations of motion and kinematic relations 
the system must obey at all times during the development of crazing 
the fundamental principle of the global conservation of energy. 

The rate of work done by external forces and all energies that enter 
or leave the material body containing crazes per unit time must 
equal the time rate of change of the internal and kinetic energies 
plus the energies associated with the formation of crazes, excluding 
chemical, electrical, etc. energies except heat. On this basis the 
global conservation of energy per unit time for a crazed visco- 
elastic medium containing u(t) number of crazes at any time t 
becomes 


W(t) = E (t) + D(t) + K (t) + Z £.(t), 

i=l (1) 

where W(t) is the mechanical work done by external load, 

E(t) is the stored elastic strain energy of the uncrazed 
portion of the medium. 
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D(t) is the dissipated energy by the uncrazed j 

portion of the medium, ! 

K(t) is the kinetic energy of the uncrazed medium. 

t 

j 

and jj 

i ± (t) is the total energy absorbed by the ith craze developed 

j 

in the medium from the beginning of craze initiation 

■i 

to current time t. jj 

In the case of relaxation it is expected that the external work j 

done W(t) approaches a constant value initially. The energy required j 

for the inception and propagation of the crazes comes from the j 

release of the stored elastic strain energy E(t) . in the case of j 

creep both the stress and strain increase and the external work I 

increases continuously which provides energy for the initiation ' 

and propagation of crazing. Each craze is an energy sink which j 

absorbs and dissipates energy from its own neighborhood. The energy | 

absorbed ? i (t) by the ith craze is characterized by the develop- j 

ment of the craze through the release of the elastic energy in the I 

surrounding medium. j 

j 

In what folio, is an attempt to reduce the problem containing j 

multi-crazes to one dealing with only one representative craze. 

Originally the whole specimen is a region subjected to an 
external stress o q , the medium is homogeneous, there is no new ! 

Phase created until crazes are developed. After the development ? 

of the new phases, subregions may be considered according to the ^ 

number of crazes developed in that region. I \ 

!j 

At any subregion it appears reasonable to consider that the j 

number of crazes developed will be proportional to the energy jj 

1 

;] 
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available in that subregion. If o (x,z,t) and e (x,z,t) are 
respectively the effective stress and the strain rate for the 
subregion, the energy available per unit volume at time t will be: 


a e (x,z,x) e e (x,z,r)dT . 


Since crazes are sources of energy sinks, several hundred times 
the amount of energy will be needed for the initiation of a single 
craze [information will be published elsewhere] as compared with 
that for the elastic deformation of the medium of tne craze size, 
essentially the craze density n(x,z,t) in this subregion will be pro- 
portional to the energy given in (2), that is 


| o e (x,z,T) e e (x,z,T)dt ~ n(x,z,t). 


In the case of creep, e e (x,z,t) is a constant after a transient 
period, then differentiation of (3) with respect to time yields: 

e a (x,z,t) ^n(x,z,t), ( 4 ) 


o e (x,z,t) 'v n(x,z,t) . 


With this important relation established, the computation of the 
time dependent craze length will be greatly simplified provided that 
n(x,z,t) is predetermined. 


IV. CONSIDERATIONS OF CRAZE PROPAGAT ION 


Assuming that each craze has basically the same configuration 

as any other craze then by referring to a central fixed (x,z) 

coordinate system, an idealized symmetrical craze in a tensile 

local effective stress field o e (t) is shown in Figure 7. By taking 

into consideration the time dependent nature of both the polymer 

medium and the craze, a craze length is designated by 2c (t) at 

time t and the time dependent stresses and displacements are 

shown in Figure 8. The half vertical distance between the two craze 

surfaces will be called the erase opening displacement and denoted 

by w(x,t) . The craze fibril domains are formed by continuous 

drawing from the unoriented bulk polymer. This drawing process 

cuases the material mass to flow from the original polymer into 

a highly oriented new phase. The stress acting on the interface may 

be called the craze envelope stress with notation o <x,t) which 

c 

is time dependent. The general state of stress in the bulk poly- 
mer is designated by o(x,z,t). 

The idealized craze structure has been considered as cylindrical 

domains of fibrils of diameter 6 f (x,t), which ma be different 

from position to position. During the crazing process the 

diameter of each fibril domain may also change with time as its 

length changes. Each fiber sustains a stress a^(x,t) which may be 

different from the craze envelope stress a (x,t) . The former 

c 

represents the true fibril domain stress while the latter is an 
average stress usually referred to as an engineering stress by 
taking the fiber domains and the voids altogether into consider- 
ation. Thus the ratio cr^(x,t)/o c (x,t) = X(x,t) equals the draw 
ratio of the fibril domains within any craze. Its inverse is the 
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volume fraction V f . The actual fibril distribution density 
function n f (x,t) is defined as the number of load bearing 
fibrils passing through a normal unit area in the craze plane at 
time t. The volume fraction of fibrils V f (x,t) inside an 
individual craze is equal to J 6^(x,t)n f (x,t) . 

The absorbed energy S(t) nucleates and creates fibril domains 
near the craze tip and transforms more material mass from the 
original phase into the new phase of oriented fibril domains, 
void formations and further deformations as well as possible 
microfractures of the highly oriented fibrils. The local 
energy balance equation per unit time for a craze at time t 
may be expressed as an energy rate balance equation with dot 
representing the time derivative: 

4<t) - f c (t) + f f (t> + D f (t) + E f <t) + K f (t), (6) 

where r c is the ener 9y required to nucleate the fibril domains 
near the craze tips, 

r f is the total surface free energy of domain fibrils 
contained in the craze, 

D f (t) is the total dissipated energy contained in the domain 
of fibrils of the craze, 

E f (t) is the stored strain energy contained in the domain 
fibrils inside the craze, 

and 

K f (t) is the kinetic energy of the domain of fibrils 
contained in the craze. 
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Aasoeiated with any craze, r f , D f , E f and can be expressed 
in terms of more elementary parameters as follows: 


I 

r c(t) 

r f (t) = 4 J 

I (x,t) y f n f (x,t)w(x,t)dx. 


c(t) 

E f (t) - 4 

J (x,t)e f (x,t)n f (x, t)w(x,t)dx. 


c ( t) 

D f (t) = 4 

6f (x,t)d f (x,t)n f (x,t)w(x,t)dx. 


c(t) 

K f (t) = 4 

1 T kf (x,t)n f (x,t)w(x,t)dx. 


where Y f is the surface free energy per unit surface area of the 
fibril domain, this energy density is assumed to be a 
material constant. 

e f , d f , and are respectively the elastic strain energy, the 
dissipative energy and the kinetic energy per unit volume of 
the fibril domain. 

Since the fibril domain nucleation rate at the craze tip must be 
proportional to the craze growth velocity c(t) , hence the energy 
required to nucleate the fibril domains near the craze tip per unit 
time becomes: 

r =4fic(t). (7) 

C C 

where B is assumed to be a constant, 
c 

The energy absorption rate of a craze |(t) can be 

obtained by performing integration of the product of the 

craze envelope stress o_(x,t) and the craze opening velocity 

c 


,il. nuif!p*rfi i_ l 


w(x,t) along the craze length, i.e. 


c(t) 


|(t) => 4 J q O q (x, t) dx 


Taking the derivatives of r f (t), E f (t), D f (t) and K f (t) 
and combining with (6), (7) and (8), one obtains the following 
equation which governs the individual craze length c(t) as a 
function of time t: 


c(t) 


” I8 f (x,t)v f (x,t)w(x,t) ] }dx 


where 


! 0 c c(t) . 


v f 2 J 6*(x,t)n(x,t) 


6 f (x,t) - — + [e f (x,t) + d. (x, t) + k-(x,t)l. 

fi f (x,t) f f (11 


Equation (9) is a nonlinear differential-integral equation of 
c(t) involving five parameters a Q , w, V f , 0 f and 0 C which form 
the basis to characterize each individual craze. This equation 
is not only describing the steady-state craze growth, but also 
governing the craze-crack transition. Once a crack is formed 
within the craze, (9) can be used to predict the crack prop- 
agation as a special case. 




Among all the parameters Involved In equation (9), moot of 
them are associated with intrinsic material proportion which are 
innensitive to the geometry of the craze and to the interaction 
of neighboring crazes. For example, at the time well before eraze- 
craek transition, to a first approximation, one may assume that 
the load-bearing fiber volume fraction V £ (x,t> and the fiber diameter 
« £ <x,t) are constants both in position x and time t during crazing 
process, it may also be assumed that the kinetic energy K f , 
dissipative energy D f and the elastic strain energy Ef are 'negligible 
Pared with the other energy term, namely the phase transformation oner 
4Y f /6 r The “ <9) reduces to the following form: 

fC(t) 

,o Mx.t.m^x, (U , 


fSjtx-t) - f o C(t) h(x,t)9w « f . ,t) d x, 


where 


and 


h(x,t) = o c {x,t) - (3 f V f , 


( 1 ?) 


5 4 6J 

which may be regarded as the energy required to rearrange a unit 
volume of the original medium into the new phase composed of oriented 
fibril domains. Equation (12) has the meaning that the craze growth 
rate c(t> is proportional to the weighted average of the craze 
opening velocity w(x,t) with weighting functxon h(x,t) as defined 
in (13). Physically (12) may be visualized as the crazs growth 
criterion derived from the first law of thermodynamics. Therefore, 
it is not surprising that (12) should be more general than the 
craze growth criteria used in [7] and 18] . 


As we have mentioned 


earlier in the introduction of this paper , Verheulpen-Heymans and 
Bauwens [7] obtained their craze growth criterion by assuming that 
the ratio of the craze opening at center (x = o) to the craze 
length c(t) be kept constant during the crazing process. This can 
be expressed in the following mathematical form in terms of our 
terminologies : 

- *1 - constant. (14) 

By taking the derivative with respect to time t, one can 
rearrange (14) in the form: 

y^t) = w (c »t) . (15) 

Comparing (15) with (12) , one can easily see that (15) is a 
special case of (12) when the weighting function h(x,t) in (12) 
is reduced to: 

h (x,t) = 6(x) , (16) 

where 6(x) is the Dirac-delta function. 

On the other hand, in order to take the strain rate behavior 
into consideration for the craze growth criterion. Argon and Salama 
[8] assumed that the craze growth rate is linearly proportional 
to the strain rate of the craze at the position a(t) which is at 
some point between the center x = o and the craze tip x = c (t) . 

This can also be expressed as follows: 

Y 2 c(t) ■ w(a(t) ,t) , 

where y 2 is a constant. 


o < a(t) < c(t) , 


(17) 


Once again (17) is seen to be a special case of (12) simply 
by taking the weighting function h(x,t) as a delta function that 

h(x,t) = 6 (a (t) ) . (18) 

Now equation (12) is a nonlinear differential-integral functional 
equation of the form [see Appendix A] ; 


f(c(t), c (t) , o Q , 8 c/ 8 f , V f , J(t), t) = 0 (19) 

• 

which is applicable for any specific applied stress o q . Among 
all the arguments, c(t) and c(t) are being sought, 0q can be 
selected, $ c , $ f , and V f must be experimentally determined, 

J(t) can be measured and t is an independent variable. Other 
qualtities such as w(x,t) must be obtained first before c(t) is 
evaluated. Some further considerations are also given below in 
connection with the analysis of this craze propagation. 


— - further conside rations of craze propagation 

In order to obtain c(x,t), w(x,t) in (12) must be obtained 
first. To determine w(x,t) the field equations governing the 
quasi-equilibrium and compatibility conditions must be solved 
together with the constitutive behavior of the viscoelastic 
media subject to proper initial and boundary restrictions. 

Using the classical correspondence principle developed in 
linear viscoelasticity and molecular orientation theory [12] the 
craze opening displacement for a single craze in viscoelastic 
medium has been presented previously [11] and the results will 
be recalled briefly here: 


w(x,t) - [C b (0)4>(x,t) + j t C b (t-T)<Mx,T)dT], 


( 20 ) 


where 


4>(x,t) = 


? r-SfS: 

x / [n 2 -x 2 ) 


( 21 ) 


g(n#t) 


■L 


° e (t) - o c (x,t) 


/(n 2 -x 2 ) 


dx , 


( 22 ) 


and 


s £-‘ 


l-s 2 v 2 

S 2 E(S) 


^-1 2 ( 1 - V 2 ) 

s 2 E(s) 


(23) 


in which the bar denotes the Laplace transform of a function 
of time, into a function in the LaPlace domain in s, i.e. 



St F(t)dt, 


/ ^ A \ 



wnere 


aenotes the inversion of the Laplace transform and v,E 
are respectively the Poisson's ratio, and the elastic modulus of 
the original unoriented bulk polymer. 1 is the natural draw 
ratio of the craze fibril domain. 

The erase envelope stress 0< .(x,t) in ( 22 ) usually depends 

upon both position and time. Based upon a number of previous 

"suits [ 9 , 11, 13 ), the stress o 0 (x,t) to a first approximation 
can be taken as follows: 


° c (x,t) * ® x (t) for o < | X | < a(t) 
o c (x,t) = a 2 (t) for a (t) < |x| < c ( t ). 


( 25 ) 


Where a(t, separates the stresses V t, and , j( t, which are both func- 
tions Of time, it ha. been shown that the use of the two-step stress 

distribution ( 9 . 11, 13, results in excellent approximation as compare, 
with actual findings. 

in addition, the introduction of step functions in the stresses 
Should reduce tremendously the mathematical complicities and computa- 
tions. This has been evidenced by the recent report on a time 
dependent theory of crazing behavior in polymers (11) . 

Because the craze envelope stress o c ,x,t) must balance the 
effective stress V t> at all times in such a way that the 
stress field within the uncrazed medium is finite everywhere 
(i.e. there are no stress singularities), then a(t) can 
be determined by the following condition!!!) . 


q c (x,t) 

(c 2 (t)-x 2 ) l ^ 2 


dx 




(t) . 


c(t) 
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By considering o^t) and o 2 (t) as follows! 


°l (t) = a i a e (t) ' 


o 2 it) m a 2 ® e (t ) • 


(27) 


where e^Cj ere constants, indeed the solution of (12) is 
greatly simplified. According to the experimental results 
in PS the stress o 2 (t) is about 10% below the applied stress 

along the majority of its length and 15% above the applied stress 
for o 2 (t) around the craze tip 114], 

Substituting (27) into (26) and solving: 


whpre 


a(t) = c (t ) cos6 , 


a . IT Vt)-0i<fc) 
9 ■ ? oJTtT=aJltT 


(28) 


TT ^~ a i 

I a 2 -a x 


(29) 


With respect to the time dependent craze density n(t) based 
upon experimental measurements at the location where craze lengths were 
monitored, the following equation was found to be quite adequate for 
representing the craze population as a function of time: 

n(t) «= A(l-e“ Bt , . (30) 

where A and B are constants. For polystyrene, referring to 
Figure 3, A was found to be 78 and B = 0.4. 

According to the linear relationship between the effective 

stress o e (t) and n(t), (see Equation (5)), the following equation can 
be established for creep: 
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a e (t) 


Ce" Bt + D 


(31) 


where c ana 0 are again constants detectable by the Initial 

tion that c e (0) = c e , and that the weighting function in (13) 

»»st be positive definite. On this basis, c ha. been found to be 

•7 N/mm o m 29.7 N/mm*, for polystyrene under the tensile 
stress c 0 - 34.4 N/mm*. This raflects the correct 

Picture that initially no erase exists and the medium is homogeneous. 

ter erases develop the erase envelope stress o c (x,t) must always 
he greater than the erase transformation energy 6fVf as shown in (13) . 
Otherwise erase will not occur or propagate. 

Apparently, information concerning the interactions among 
neighboring erases can be effectively taken care of through the 
application of the effective stress , e (t) as a function of the 


rate of change of the local erase number. 
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VI . RESULTS 


Using computer aided numerical methods (A15) a reduced form of 
(Al) or (19) or (12) has been solved. For polystyrene, only 
three constants are now needed. They are 
3 C = 0.7 x 10" 6 J/mm 2 

0£ = 0.112 J/mm* 

and V f = 1/4 


These values are obtained on the basis that 6 f = 15 nm, Y 
J/mm 2 and X = 4 which have been experimentally found [15] . 


f = 0. 42» 10“ 6 
Note that 


V f 5 1/A - Using these together with the creep compliance J(t) as 
determined m (B2) Appendix B, the following results have been obtained 
The craze envelope stress distribution c c (x,t) for 
polystyrene is given in Figure 9 for a specific craze of length 
just below 150 ym. Based upon various considerations, as 
stated before, a fairly accurate craze propagation behavior 
should prevail if all the subregional crazes are consistently 
stressed, indeed this is the case, the computed time dependent 
craze length is given in Figure 10 in which one normalized curve 
is shown to represent the typical craze length-time behavior. Com- 
pared with the experimental data as shown in Figure 6, a fairly good 
agreement with the calculated results is obtained. 


It nay be quite significant that the tine dependent single 
erase length in a viscoelastic medium has been found to propagate 
to infinity in a finite time if only one erase is considered 111) 


According to the current analysis, when craze interactions a 
taken into consideration, however, the craze length as a 
function of time propagate to a finite value as time goes to 
infinity. 



VII * DISCUSS ION 


(1) Craze Propagation 

The craze propagation computed analytically agrees fairly 
well with the experimental result. (Figure 10) . it was found that 
the craze growth is decelerated with time. This means that the 
craze length increases with time and its velocity decreases 
with time. From Eguation (12) , it is observed that craze arrest 
(i - 0) occurs when either a c - 6fVf or 9w(x,t)/9t - 0. 

The former occur, when the energy absorption rate by crazing is 

used up for drawing new fibrils out of the matrix, while the 

latter is possible only if molecular entanglements enhance the stiffnes, 

of the fibril domain such that further opening becomes difficult. 

(2) Craze Number 

As shown in Figure 11 it was found that the number and 
length of crazes developed in creep are strongly dependent on the 
magnitude of the external load. At high stresses a large number 
of small crazes become visible after a relatively short time. 

while under small stresses a small number of large crazes become 
visible after a long loading time. 


(3) Local Effective Stress and Rate of Craze Number 

It is interesting to note that the local effective stress is 
linearly associated with the rat. of change of the local craze 

number. This relationship seems to be quite useful for analyzing 
craze interactions. 
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APPENDIX A 


On Craze Length e(t) 

By substituting (25) and (27) into (12) one obtains: 

a(t) 

la l° e tt> • B f V f ) { 0 - W - ( >4^ dx 


+ la 2 a e (t) - 0 f V f ] 


c (t) 


3w(x, t) 


Tt 


dx 


3 c c(t) 


a (t) 


(Al) 


Further, by substituting (25)and(27) into (21)and(22), <Mx,t) is 
found as follows: 


4>(x,t) = i(o 2 -a 1 )o e (t) [(x+ycttWcosh-'SiJ^llx 

x+ c(t) 


- (x-yc(t) Jcosh” 1 


(A2) 


where 


Y * cose, 


<Mx,t) 

■x=a(t) 


(a 2 -o 1 )o e (t)c(t)I 2 . 


(A3) 


with 


° - cose In sece. 


(A4) 


1 


o 


<P(x,t)dx = (« 2 -« 1 )a e (t)c 2 (t)l 1 , 


(A5) 


re 


I x * cos*d( | tane - i tane + l in seed), 
t c(t) 

Ja(t) 4,(X,t)dX = (o 2-«i)^^ c2 (t)I 3 , 


(A6) 


(A7) 


e 


Z 3 * ? cos 2 0 (0tan0 - 2 In sec0) . 


(A8) 


d f a < fc > 

J Q <Hx,t)dx = (a 2 - a;L ) [d e (t)c 2 (t) + 2a 


dt L(t) <,>(X,t,dX * (a 2" a l ) ta (t)c 2 (t) + 2a 


c(t) 


e (t)c(t)e(t)]z 
(A9) 


a(t) 


e (t)c(t)c(t) ]l 
(A10) 


Using the above results anri tons - 

esuj.es and (20 ) , one gets the following; 

f a ( t )g W ( x# t) ! 

J 0 3t ** 58 T=l HC^OJI^tJcMt) 


+ C b<°H2l 1 -i 2 )a e (t)c(t)c(t) + C b (0)l 1 a e (t)c 2 (t)). 
(a 2 "°i) + R x (t)} f 



where 


fc(t) 

W)^?' ,t)dx “ r=r nc b (o)i 3 a 0 (t)c a (t) 


+ c b (0) (2I 3 +I 2 )a e (t)c(t)i(t) 
+ R 2 (t)}, 


+ < = b (0)I 3 o e (t)c 2 (t)] (a 2 - ai , 

(A12) 


; a (t) f t m 

R l(t) = I C*(t-T)«Mx,T)dTdx, (A13) 

J o * o 

f c(t) ft„ 

R 2 (t) “I C b (t-TH<X,T)dTdx, ( A14) 

J a(t) J o ° 


and X is the fibril draw ratio. Substituting (All) and (A12) into 
(Al ) , finally one obtains: 


c(t) 


K 1 (t) c 2 (t) + K 2 (t) 

B c - K 3 (t)cTt) 


(A15) 


where 


K l (t) " X=T (a 2 “ a l ) tC b (0)5 e (t) + C b (0)a e (t)] . 


lI l L l (t ) + I 3 L 2 (t)1 (A16) 

K 2 (t) “ X=r lL l (t)R l (t) + I*2 (t)R 2 (t)] (A17) 

K 3 (t) " T=r< a 2“°l )C b (0)I(2l l“ I 2 )L l (t) + < 2l 3 +I 2 )L 2 (t,la e (t) ' 
with < A18 > 

L l (t) " a l a e (t) “ 8 f V f' (A19) 

L 2 (t) = o 2 0 e (t) “ 6 f v f (A20) 
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* APPENDIX B 

On Creep Compliance J(t) 

Before any calculation can be accomplished, the constitutive 
behavior of the medium in creep must be obtained and properly 
represented for computation. 

Taking the experimental polystyrene sheet as an example, 
Figure 12 shows the experimentally measured creep compliance J(t). 
The quantity C^(t) is closely related to this creep compliance. 

By assuming for polystyrene that the Poisson's ratio v = 0.395, 
a constant, then 


C. (t) 1 2 * 1 ~ v>2 - ) ■ 2(l-v 2 )e?f J(s) 

° s 2 E(s) 


2 (1-v 2 ) J (t) . 

(Bl) 


For any arbitrary times, a very good approximation may be obtained 
using a linear viscoelastic spectrum form [16] : 

m /T • 

J(t) * J n + E J. (1-e w x ) 

0 i=l 1 (B2) 

where J^'s andt^'s are respectively constant material compliances 
and retardation times. With the measured creep compliance curve 
available, these quantities can be obtained using a curve fitting 
technique. 

For polystyrene under the creep stress a Q = 34.4 N/mm , the 
following data are found to represent a close fit: 


4 
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m 


m - 4 j^ a 

J Q = 0.307 x 10" 9 mm ? /N 
a 0.071 x 10“° mm ? /N 
J 2 = 0.062 x IQ’ 9 nun 2 /N 
J 3 = 0.045 x 10- 8 mm 2 /N 
J 4 = 0.031 x IQ" 9 mm 2 /N 


t x ™ l.Q/h, 
t 2 a 10.0/h, 
t 3 <=> 80.0/h, 
t 4 o 110. 0/h, 
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